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INTRODUCTION 

The turbulent boundary layer (TBL) pressure field is an important source of cabin noise during cruise 
of high subsonic and supersonic commercial aircraft [1-3, 6], The broadband character of this excitation 
field results in an interior noise spectrum that dominates the overall sound pressure level (SPL) and 
speech interference metrics in the forward and midcabins of many aircraft. In the authors' previous 
study [9], sound transmission through an aircraft fuselage, modeled by two concentric cylindrical 
sandwich shells and excited by a TBL statistical model was investigated analytically. An assessment 
of point and global structural vibration levels and resulting interior noise levels was obtained for 
different TBL models, flight conditions and fuselage structural designs. However, due to the 
complication of the shell structure, the important noise transmission mechanisms were difficult to 
discern. Previous experience has demonstrated that a fundamental understanding of the range of modes 
(or wavenumbers) generated by the TBL source both in the structure and the acoustic cavity is key to the 
development of both active and passive control technologies. 

In an initial effort to provide this insight, the objective of this paper is to develop an analytical 
model of sound transmission through a simple unstiffened cylindrical aluminum shell excited by a TBL 
pressure field. The description of the turbulent pressure field is based on the Corcos formulation [2] for 
the cross-spectral density (CSD) of the pressure fluctuations. The coupled shell and interior and 
exterior acoustic equations are solved for the structural displacement and the interior acoustic response 
using a Galerkin approach to obtain analytical solutions. Specifically, this study compares the real 
part of the normalized CSD of the TBL excitation field, the structural displacement and the interior 
acoustic field. Further the modal compositions of the structural and cavity response are examined and 
some inference of the dominant mechanism of noise transmission is made. 

MATHEMATICAL MODEL 

A schematic of the system configuration considered is shown in Fig. 1, a finite cylindrical shell 
having hard end caps with radius R, length L, and wall thickness h immersed in fluid media. The mass 
density and sound speed are denoted as pi and c\ and p 2 and C 2 for the external and internal fluids, 
respectively. Turbulent boundary layer pressure fluctuations convected by the flow in the external fluid 
along the streamwise direction excite the surface of the shell. A detailed analysis of the pressure 
fluctuations models has been presented in the references [2, 3, 6] and will not be repeated here. W e 
shall instead, simply quote an expression for the spatial-frequency CSD of the pressure fluctuations S p 
as following, 

Sp(£, i), co) = <P(co) A(%,co) B(ri,co) e im t' Uc 


( 1 ) 
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where i = V-l, ©the frequency, <; = xi-x 2 the streamwise separation along the longitudinal x direction 
and r) = R(9 i-9 2 ) the crossflow separation along the circumferential 9 direction, U c the convective 
velocity of the flow, and <P the power spectrum. Quantities A and B are exponential decaying functions 
along g and ij. 

Taking a Fourier transform of a Helmholz equations for the external and internal fluids and 
Donnell-Mushtari shell equation yields, 

(V 2 + icofh + k\) pi = 0, (V' + i©/L + ki) pi = 0 (2) 

L(©)[w] = pi-pi-p ^ (3) 

where V 2 = d(r d/dr)dr + ( d 2 /G9 2 )/r 1 + d 2 /dx 2 in which r represents the radial cylindrical coordinate; 
k\ = co/c i and k 2 = co/c 2 are wave numbers, pi and pi acoustic damping, and, pi and pi acoustic pressures of 

the external and internal fluids, respectively; p is the pressure generated by the TBL, L is a 
differential operator whose explicit form can be found in text book [4], and zv is the radial displacement 
component of the shell. 

The flexible boundary conditions applied at both the outer and inner surfaces of the shell relate the 
pressures pi and pi to the shell displacement zv 

= pi© 2 w, = p 2 a> 2 w (4) 

dr l=* dr r=R 

A simply supported shell having hard end caps is represented by the following boundary 
conditions. 

u, - d w - n ^p2 = n at r = n r 


d 2 w _ 0 dp 2 

dx 2 dx 


0 , at x = 0, L 


Following the previous analysis [9], the displacement zv is expanded in the orthogonal base 
functions X mn 

oo oo 

W = WmnX m „(x,0 ) - Wmn sm(m7rx/L) cos(nG) (6) 

m = l, n = 0 m = 1, n = 0 

where the modal coefficients w mn are unknown. Similarly, expand the turbulent pressure loading p as 


p = y Pmn X mn (x,9 ) , p m „ = [I p( x ,9,(o)X mn (x,0)dxd0] 

m = 1, n = 0 7lL - 0 JO 

where £ n are Neumann factors, £ n = 1 for n= 0 and £ n = 2 otherwise. 

It is convenient to expand the pressure pi and p 2 which satisfy Eq. (2) as 


Pl — Plmn Xmn ■> P2 — P2mn J n(&2ml r ) Xmn (8) 

m = l, n = 0 m = 1 , n = 0 

2 

where Hk is a Hankel function of the second kind and /„ is a Bessel function of the first kind of order n, 
respectively, a 2 m = k\ - (mK I L) 2 - icopi and a 2m = k 2 - {mn /L) 2 - iojpi, and pi mn and pimn are unknown 
modal coefficients. 

Substituting Eqs. (6) to (8) into (2) through (4) and then solving these coupled equations 
simultaneously for w mn , pimn, pimn yield, 

{ W mn-> pimn ? P2mn } - {H wmn-f H p i mn 5 Hp2mn} Pmn ( 9 ) 

where 

{H wmn •> H p i mn 9 Hp2mn} — { l/(i Of), — Zmn! J mn, ~ Zmn / Hmn } / (Zmn +Zmn ~ Zmn) (10) 

Xmn — ip 1 lOH m ,J(C\,n H >t!n) . Xmn — ip 2 (OJmndOlm J mil) (H) 

The primes in above equation denote derivative with respect to the argument and X,„ n is the shell 
modal impedance. 

The cross-spectral density for the wall displacement response and interior pressure takes the form 
S w - y \H wm „\ 2 sm(mnxi/L) sinimnxi/L) cos(h0i) cos(n9 2 )S Pmn (12) 

m — 1 , n = 0 


Sp2 — \Hp2mn\ \Jn(&2mT’l) Jni.^2m^2) COS(wt7Tri/E) COS(wt7Df2/E) COS(/10i) COS (jl02)Sp mn (13) 

m - 1, n = 0 
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where * denotes the complex conjugate, 0 <n, >~2 <R, 0<xi, xi <L, -k< 9\, Qi<n, and the modal spectral 
density of the generalized forces S Pmn is defined by 

S Pmn = Re I 4<P(oo) £„ (7iLR\Y l (a/T 1 lx[l-(-l)V^ x { 1+ [n /f yS/? i )] 2 } 1 x 

{ \+[tnnl(a'L)] 2 +2[\-(-\) m e- a L ]x(mn)\aLy i }x{ \+[trml(aL)f}~ 2 ' ' 

In Corcos model, a ' = (0.11+ i)w/lf c and ft - 0 .77o}/U c . 

It should be pointed out that the modal basis function for the pressure pi, Eq. (8), does not satisfy 
the force boundary conditions in Eq. (5). Therefore, it is necessary to expand the orthogonal base 

functions X mn in terms of cos(mKxIL) cos (n6). The result gives the modal coefficient H p i mn as following 


Hplmn 


i 

j = o,j*m n (m -j ) 


(15) 


Since both the excitation and cylindrical shell are symmetric in 9, integration of Eqs. (12) and (13) 
yields the average CSD 


{Sw, Sp 2 } — Sp2mn\ ~ { \H wmn \ , \H p 2 m „\ }Sp m J 2 

m = l, n = 0 m = 1 , n = 0 

Summation of the components either m or n in above equation yields the modal CSD 


{Sw m , S p 2 ,„} — Sp 2mn } , {5w„, S p 2„} — S p2mn } 

n = 0 m = 1 


(16) 


(17) 


RESULT AND DISCUSSION 


The geometry and material properties of the shell chosen herein are taken from reference [5] for 
this study. The aluminum shell has the overall dimensions of R - 0.83m, L = 3.65m, and h = 0.0017m. 
Young's modulus and Poisson's ratio are E = 72GPa and ll = 0.3, respectively and the mass density is p = 
2740kg/m 3 . A structural loss factor is modeled as a small imaginary part of both Young's modulus 
£(l+i0.05) and a corresponding shear modulus. The value of acoustic damping is chosen as 0.005 
providing representative damping levels. It is assumed that the aircraft is in flight at 35,000ft 
altitude (pi = 0.3795kg/m 2 , c\ = 296.556m/s) and is pressurized at 10,000ft (pi = 0.9041kg/m 3 , ci = 
328.558m/s). The free stream velocity of the turbulent flow is taken from Maestrello [6] for LG = 
425.48m/s and then the convection velocity is defined by U, = 0.8LG. 

The result presented in the following will be two parts: Part one will address the coherence 
functions of the two-point spatial correlation of the excitation, shell response (displacement), and the 
cavity response (pressure). It is pointed out [8] that the modal spectra for the shell and cavity 
responses are more useful from the viewpoint of developing noise control systems than the highly 
detailed physical variables. Therefore, part two will examine the modal response of the cross- 
spectral density for both the shell and acoustic cavity responses. A dominant (or highly excited) mode 
is indicated by a dark square. The values in all the following plots are normalized by its own maximum 
and only the real part is considered. 

Figures 2 and 3 show the CSD of the excitation (TBL pressure fluctuations, S p ) and the corresponding 
shell and acoustic responses ( S w and .S' ; , 2 ) versus the spatial separation in the streamwise q/L and cross- 
streamwise q/(27rR) for co - 300Hz and 2500 Hz, respectively. The reference coordinate is selected as x 2 
= 0.48L to avoid the nodes at which some modes can not be excited. Therefore, a slightly asymmetric 
behavior is shown in Fig. 2. Since the excitation is independent of the circumferential coordinate 
selected, the reference coordinate is 02 = 0. Inspection of Fig. 2 reveals that the envelope of the shell 
response decays with an increase of the spatial separation, in a similar manner as does the excitation. 
The rate of decrease depends on the frequency of the forcing function. The results illustrate that the 
shell response at two widely separated points is generally correlated at low frequency a> - 300Hz and is 
uncorrelated at high frequency a> = 2500Hz. However, for the acoustic response, we do not observe a 
similar decaying relationship as for the shell response. The new feature can be explained by the 
following: When the TBL pressure fluctuations excite the shell at a particular frequency, many shell 
modes are excited as resonant modes. However, only some of these shell resonant modes are efficiently 
coupled to those cavity modes near resonance. Therefore, the dissipation of the acoustic CSD along the 
spatial separation can not be seen when an insufficient number of cavity modes are excited. If the shell 
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radius is increased (actually increasing the excited cavity modes), a stronger decay of the CSD is 
observed. 

Figure 4 shows the modal cross-spectral density of the shell response S Wn and the acoustic response 
S P2n versus the circumferential mode n and frequency co. High mode levels occur in the shell response in 
a line between 13<n<22 over the range of frequencies plotted, as shown in Fig. 4(a). A strong filtering 
effect is observed in comparing the cavity mode response and the shell wavenumber mode response. 
Most low order acoustic cavity modes between 1 <n<7 are strongly excited, but most of the high order 
circumference modes are reduced in the cavity response. 

Figures 5 and 6 display the modal cross-spectral density of the shell and acoustic responses S wmn and 
S P 2 mn for co - 100Hz and 300Hz, respectively. There is strong shell response over a range of resonant 
shell modes and this range expands to higher order modes with an increase of frequency. This 
demonstrates that with increasing frequency, an expanded locus of resonant modes is responding. For 
the acoustic response, these resonant shell modes drive a forced response of the corresponding acoustic 
modes and in addition excite some low order cavity modes into resonance. Similar behavior was 
observed for the modal structural and acoustic response from point force excitation [8]. This shows that 
the CSD of the cavity response is not only characterized by the "forced" modes generated by the shell 
vibration, but also by the cavity modes, especially for high frequency. 

CONCLUSION 

In this paper, an analytical model for prediction of the CSD of both the shell and acoustic 
responses of a finite cylindrical shell due to TBL excitation based on the Corcos TBL model was 
developed. The real part of the normalized CSD at two-points versus a spatial separation of the 
forcing field and response fields are calculated and presented. Modal behavior is examined in the low 
frequency range (100 Hz to 300Hz). The following conclusions can be drawn: 

The CSD of the shell response are generally correlated at low frequency and the envelope of the 
correlation function decays with an increase of the spatial separation. The decay occurs much more 
rapidly at higher frequencies. However, the correlation function of the CSD of the cavity response 
does not decay with an increase of the spatial separation and frequency. With wavenumber filtering, 
highly excited cavity modes are reduced in number compared with the shell modes. Most of the cavity 
mode response is forced by the shell modes. However, some low order cavity resonant modes are 
excited. 
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Figure 1. Schematic of system configuration. 




.. The normalized cross-spectral density versus the spatial differences i;/L in the streamwise direction: 

(a) a = 300Hz; (b) co = 2500Hz. 




Figure 3. The normalized cross-spectral density versus the spatial differences £/L in the cross-streamwise direction: 

(a) co = 300Hz; (b) co = 2500Hz. 
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Figure 4. Modal cross-spectral density versus circumferential mode and frequency: varying from 

0 to -30dB. (a) shell response S w „; (b) cavity response S P2n . 
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Figure 5. Modal cross-spectral density S Wmn versus longitudinal and circumferential modes m and n: 

varying from 0 to -30dB. (a) co = 100Hz; (b) <o = 300Hz. 
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Figure 6. Modal cross-spectral density S P2m „ 
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Longitudinal modal number, m 
longitudinal and circumferential modes m and n: 


varying from 0 to -30dB. (a) co = 100Hz; (b) co = 300Hz. 
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